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In this paper, we present how to use an interval arithmetics framework based on free algebra 
construction, in order to build better defined inclusion function for interval semi-group and for its 
associated vector space. One introduces the ?/)-algorithm, which performs set inversion of functions 
and exhibits some numerical examples developped with the python programming langage. 
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I. INTRODUCTION 

Many problems dealing with real numbers can be solved with intervals computationd^'^'^MlSI^ such as set inversion, 
functions optimization, linear algebra. Their key stone is the construction of the inclusion function from the natural 
one. Some approaches have developped methods to circumvent it with using boolean inclusion tests, series or limited 
expansions of the natural function where the derivatives are computed at a certain point of the intervals. Nevertheless, 
those transfers from real functions to functions defined on intervals are not systematic and not given by a formal 
process. This yields to the fact that the inclusion function definition has to be adapted to each problem with the 
risk to miss the primitive scope. 

We propose here a clear and simple scheme to build inclusion functions from the natural one. It is based on the 
arithmetic developped by one of the author in his PhD thesid^, called in this paper free-algebra based arithmetic. 
This article reminds first the definition and characteristics of the intervals semi-group IM and the construction of 
its associated vector space IM. We explain how to get an associative and distributive arithmetic of intervals by 
embedding the vector space into a free algebra. After that, one shows from an example how to generalize the 
inclusion function construction for the semi-group and the vector space. This permits to present the ^/'-algorithm 
in order to perform set inversion which is based on this construction and probability calculations. We ends with 
numerical applications examples developped in pythort^^ programming langage for set inversion. 



II. THE BANACH SPACE IK AND THE ASSOCIATIVE ALGEBRA Aa 

An interval X is defined as a non-empty and connected set of real numbers. It can be represented as a point in the 
half-plane of K^, Vi — {{a,b) S E^, a < b}. The set V2 = {(1, &) G M^, a > b} is the set of anti-intervals as shown on 
figure [I] One writes real numbers as intervals with same bounds, Va e M , a = [a, a]. We denote by IM = Vi the set 
of intervals of M. The semantical arithmetic operations on intervals, called Minkwoski or classical operations, 
are defined such that the result of the corresponding operation on elements belonging to operand intervals belongs to 
the resulting interval. That is, if o denotes one of the usual operations — , *, /, we have, if X and Y are bounded 
intervals of M, 

XoY ^{xoy / xeX, yeV}, 
But IM is only a semi-group for the addition and no substraction can be defined for it. 



A. Modal interval arithmetic 

An algebraic extension of the classical interval arithmetic, called generalized interval arithmeticP^ has been pro- 
posed first by M. WarmuJ^. It has been followed in the seventies by H. -J. Ortolf and E. Kaucheil^^. In this 
former interval arithmetic, the intervals form a group with respect to addition and a complete lattice with re- 
spect to inclusion. In order to adapt it to semantic problems, Gardenes et al. have developed an approach called 
modal interval arithmetictli^HII!. g. Markov and others investigate the relation between generalized intervals opera- 
tions and Minkowski operations for classic intervals and propose the so-called directed interval arithmetic, in which 
Kaucher's gene ralized intervals can be viewed as classic intervals plus direction, hence the name directed interval 
arithmeticP^E^. In this arithmetic framework, proper and improper intervals are considered as intervals with sigrP^. 
Interesting relations and developements for proper and improper intervals arithmetic and for applications can be 
found in litteraturtpSHU) 

Our approach^^ ^% that we remind below in this artice, is similar to the previous ones in the sense that intervals 
are extended to generalized intervals. However we use a construction which is more canonical and based on the 
semi-group completion into a a group, which permits then to build the associated real vector space, to get an analogy 
with directed interval arithmetic and to perform differential calculus and linear algebra naturally. 
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B. New interval arithmetic 



Let IM be the IM completed group. An element of IM is written {A, B) where A, B e IM and the equivalence 
relation being {A,B) - {A',B') if and only ii A + B' = B + A'. For example ([3, 4], [1, 4]) - ([2, 0], [0, 0]) = ([2,0],0). 
Any equivalence class can be written as (AT, 0) or (0, X) or (a, a) with A G IM or a e M. We denote X\y,X, y £ M, 
the difference into the abelian group IM. The external multiplication 



X I! 



given by 



a ■ {X, 0) = (a • X, 0) if a > 
= (0,-a- A) if a < 



(1) 



(2) 
(3) 
(4) 



and 



a - (0,A) = (0,a- A) if a > 



= (-a • A,0) if a < 



(5) 
(6) 
(7) 



provides IM with a 2-dimensional M vector space structure. A basis is given by the family {([1, 1], 0), ([0, 1], 0)}. The 
following map 

||.||:m — ^M+ (8) 

(9) 
(10) 

(11) 



(A,0) ^ ;(A) + |c(X)| 
WX) ^ l{X) + \c{X)\ 



with respectively 1{X) and c(A) the width and the center of the interval A, is a norm. It is easy to prov^^^'^ that IM 
is isomorphic to M^ which is complete, this yields to the fact that this norm endows IM with a Banach space structure. 
Moreover, one can define an order relation < on IM which is compatible with this norm : A, F e IM, A < F if 



A r and c(A) < c(A) 



(12) 



or 



A C r and 1{X) < 1{Y) 



(13) 



Due to the fact that every intervals can be compared to each others, it yields that this partial order relation can be 
extended to a total one. 

In order to define a natural multiplication on IM, we consider the following 4-dimensional associative algebra A4 with 
basis : 



and product 



' ei 


= (1,1), 




= (0,1), 


63 


= (-1,0) 
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= (-1,- 
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(14) 
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Figure 1. Intervals and anti- intervals can be respectivelly represented as points in half planes Vi and Vi- This implies naturally 
the representation of interval vectors from IR. 



Those vectors represent the intervals [1, 1], [0, 1], [— 1, 0] and [—1,-1] (figure [T]). The product table is obtained 
simply with the mean of Minkoswki product. For example 62-63 = [0,1] • [—1,0] = [min(0, — 1), max(0, — 1)] = 
[— 1, 0] = 63. One can define the sectors Vij, with j — 1, 2, 3 such as 

j = 1 : (a, b) = aei + [b — 0)62, if a > 0. (15) 
j = 2 : (a, b) ^ -063 + 6e2, if a < and 6 > 0. (16) 
j = 3 : (a, b) = -hci + {b- a)e3, if 6 < 0. (17) 



Thus any interval X (i.e any {X, 0) G IK) is represented by positive components vector of A^. Let's denote this 
embedding (p : M — > A4. The representation of (0, X) is then given by (p{{0, X)) — ~(p{{X, 0)). Let F be the linear 
subspace of A4 spanned by 61 — 62 + 63 and 61 + 64. Any vector J^i ^a^i € A^ is equivalent in A^/F to a vector 
with positive components if 012 + 03 > 0. In this case, such a vector corresponds to an interval of IM and one has the 
bijection : IM — > V where V is the convex part of Aa/F corresponding to 0:2 + as > 0. This defines a distributive 
product in IM by 

X-Y = w\v'{X)-p'{Y)) (18) 

where Lp' = tpoip and tp is the canonical embedding from Ai to A4/F as shown on Fig. ([2]). Division between intervals 
can also be defined with solving X ■ Y = {1,0, 0, 0) in A4 or in an isomorphic algebra. 

One can remark that : 

• Distributivity is satisfied for multiplication and division according addition and substraction 

• Usually, one has X oY C X ■ Y where X oY is the classical Minkowski product. More precisely, within this 
algebra, if X and Y are not in the same sector Vij, or are in the same sector Vn or Vis, one gets the Minkowski 
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product. If they are belonging both to 7^12, the result is greater than the Minkowski one. If one needs a result 
X ■ Y closer to X o F or even equal to, one can embed IM in higher dimension algebras, such as A^^A^.A^, 
and so on. This former one is sufficient to recover the distributivity and equality to Minkowski productP^. 

• It can be naturally extended to M. 

• The product satisfies the monotony property. 

• With the canonical injection X M- {X, 0) from IM to IM, one can consider that IM C IM. 




Figure 2. Interval embedding in the algebra Aa,. 



III. INCLUSION FUNCTIONS 



It is necessary for some problems to extend the definition of a function defined for real numbers / : M — > M to 
function defined for intervals [/] : IM — >■ IM such as [/]([a, a]) — f{a) for any a € M. It will be convenient to have 
the same formal expression for / and [/]. Usually the lack of distributivity in Minkowski arithmetics doesn't give 
the possibility to get the same formal expressions. But with the arithmetic we have presented, there is no data 
dependency anymore and one can define easily inclusion functions from the natural one. For example, let's extend to 
intervals the real functions fi{x) = x'^ — 2x + 1, f2{x) = {x — 1)^, faix) — x{x — 2) + 1. Usually, with the Minkwoski 
operations, the three expressions of this same function for the interval X = [3, 4] are [f]i{X) = [2, 11], [/]2(-'^) = [4, 9] 
and [f]^{X) = [6, 12]. Data dependancy occurs when the variable appears more than once in the function expression. 
The deep reason of that is the lack of distrib utivity in Minkowski arithmetics. But within the arithmetic developped 
in Ai or higher dimension free algebra^^HHl^ t^jg problem vanishes. For example : with X = [3, 4] and since X S T^n, 



ip{X) = (3, 4 - 3, 0, 0) = (3, 1, 0, 0) 
Since ei = 1, and with means of product table (14) one has 



3ei 



62- 



(19) 



^([/]i(X)) = (3e 



2(3ei + 62) + 1 
9e? + 2 • 3eie2 + - 2 • 3ei - 2e2 + 1 
9ei + 6e2 + 62 — 6ei — 2e2 + Ci = 4ei + 5e2 
^([4,9]), 



(20) 



^{[fUx)) 



(3ei 
(2ei 
46? H 
4ei H 
4ei H 
^([4,9] 



^62- 1)" 
4eie2 + 63 

4e2 + 62 
562 



(21) 
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and 



vilfUX)) = (3ei + 62) • (3ei + 63 - 2) + 1 

= 9ei + 3eie2 — 6ei + 3ele2 + 62 — 2e2 + Ci 
= 4ei + 3e2 + 3e2 + 62 - 2e2 
= 4ei + 5e2 

= ^([4,9]). (22) 

Thus, = [/]2(X) = [fUX) = [4,9]. 

On the other hand, the construction of the inclusion function depends on the type of problem one deals with. If 
one aims to perform set inversion for example, it has to be done in the semi-group IM due to semantic definition of 
substraction. But, the substraction is not defined in IM. This problem can be circumvented by replacing it with an 
addition and a multiplication with the interval 64 = [—1,-1]. This maintains the associativity and distributivity of 
arithmetics and permits to introduce a semantical substraction. For example : if f{x) = x'^ — x = x{x — 1) for real 
numbers, one defines [/](-'f) = + €4^ ■ X. One reminds the product [—1, —1] • [a,b] is equal to [—b, —a]. Due to 
the fact that the arithmetics is now associative and distributive, one doesn't have data dependancy anymore and 
= X'^ +64- X = X ■ {X + 64). The last term corresponds to the transfer of x{x— 1). Division can be transfered 
to the semi-group in the same way by replacing ^ = x^^ with X'^'^. 

Taylor polynomial expansions, differential calculus and linear algebra operations are defined only in a vector space. 
Therefore the transfer for the vector space is done directly. This permits to get infinitesimal intervals with the 
substraction and to compute derivatives. This is of course not allowed and not possible into the semi-group. From 
IM to the vector space IM, f : x —x is transfered to [f] : X = {X,0) 1— )■ 0) = \X. This means that 
[a, b] substraction is the anti-interval [—a, —b] addition. One of the most important consequence is that it is pos- 
sible to transf er s ome functions direc tly to the intervals. For example, it is easy to prove analytically in IM that 
[exp](([a, 6], 0)) = ([exp(a), exp(6)], 0) with means of Taylor expansion. 



IV. PROBABILIST SET INVERSION : i/'-ALGORITHM 
A. Flowchart 

One presents an efficient set inversion method whose flowchart is very simple. One of the powerful application 
of interval calculus is the set inversion of a real-valued function defined on real numbers. It finds a huge number 
of applications in automation, optimization, resolution of algebraic and differential equations. The mathematical 
modelization of this problem is the following as shown on Fig.([3| : let's note / : M" 1— >■ M^ a function for a physical 
system, which is required to be surjective only, TZ C M" the set of adjustemcnts, and V CM.^ the set of performance 
of a system. Set inversion consists of the computation of 5 = f^^{V)r\Tl. For semantical reasons, one has to perform 
this set inversion within the semi-group IM. 

Some interesting and powerful methods using intervals have been developped those last years, especially SIVIAE21j 
Set Inversion Via Interval Analysis. But the inclusion function being not well defined in the semi-group with the 
Minkowski arithmetic, SIVIA uses boolean inclusion tests and finds accepted, rejected and "uncertain" domains. 
With the algebraic arithmetic, one doesn't need boolean tests since the inclusion functions are well-defined. Thus, 
we propose the "^-algorithm (Probabilist Set Inversion) inspired from SIVIA but without boolean tests and with a 
conditional probability calculation and domain bissections. This yields to accepted or rejected domains only. We are 
interested to compute the following conditional probability 

p{x) = p{[f]{x) c V I f{x) e [f]{X), Vx e x) 

^ mes{[f]{X)C^V) 
mes{[f]{X)) 
mes{y n V) mes{X) 
mes{y) mes(y) 



(23) 
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Figure 3. Set inversion modelization with intervals 



PSI of f(x,y)=(x'^2+y^2,x+y) for performance set [1,2]X[1,4] 
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Figure 4. i/;-algorithm for /i 



where mes is the Lebesgue measure in W (length, surface, ...). If this probabihty equals 1 then the set is added to 
the list of solutions. If it is zero the set is rejected and removed from the list of interval candidates. If the probability 
is such as p{X) s]0, 1[, then X is bissected and ■(/'-algorithm applies the same procedure recursively for the resulting 
intervals until the size is lower than a fixed size resolution of the intervals or until the sets are accepted or rejected. 
Since ^/^-algorithm creates sequences of decreasing intervals which are compact sets, it is obvious that ^/^-algorithm 
converges to fixed points probabilities which are simply and 1. In fact, we consider a sequence of compact sets 
{Kn}nm satisfying d{Kn) < d{Kn+i) where d is the diameter of the compact set. If Kn n Km = 0, for any n ^ m, 
the sequence is convergent and the limit is the empty set. Then, one has just to consider the sequence of bissected 
sets. 



B. Numerical applications 



We have developped a numerical library for python environnem entP^ called yet intervalJib. It is a pure numerical 
implementation performing the basic arithmetics presented abovd24|25j ^j^jg library aimes to give simple and opti- 
mized routines to perform interval calculations based on the algebraic arithmetics. One gives in this section some 
numerical application examples of -i/j-algorithm in order to illustrate how it can treat usual inversion problems and 
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Figure 5. t/j-algorithm for /2 



PSI Of f(x,y)=(x'"2-y'"2,y/(l+x)) for performance set [0,5]X[-4,4] 
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Figure 6. t/j-algorithm for /s 

build well-defined inclusion functions. 

Let's define the non-linear functions fi : with respectively adjustements and performances sets TZi, Vi : 

h{x,v) = {x-' + y^x + v), 7^l - [-1,2]2, Vi = [1,2] x [1,4] 
h{x,y) = {x^ -V\x-y), 7^2 = [-3,3]^ T'a = [0, 4] x [-2, 2] 

h{x,y) = {x^ - y\ ^), = [0,6] x [0,10], P3 = [0,5] x [-4,4]. 

(24) 
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PSI for f(x,y)=(x'^2-y'^2*exp(x)+x*exp(y),x*(x+y)-y'^2) 
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Figure 7. t/j-algorithm for /4 
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Figure 8. Accepted and candidate intervals total ratios for numbers and surfaces according iterations in ?/)-algorithm for fi 
with 10~* spatial resolution. 



Those examples have been chosen to give examples of addition, substraction, product and division transfers from 
M to IE, and to exhibit the difference between the usual Minkowski arithmetic and the algebraic on^^lF [HI, The 
calculations with the 7/;-algorithm are shown on figures |4j |5j [6] and [7| The convergence to or 1 probabilities only, 
shows that inclusion functions are well constructed and that the interval arithmetic is robust. The last example 



fi{x,y) = {x^ - 2/^ ■ exp(x) + X ■ exp{y),x ■ {x + y) ~ y^), 7^4 = 7^4 = [-5,5^ 



(25) 
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presented on figure [7] shows clearly that the ^/'-algorithm implemented in the algebraic arithmetic we use is not data 
dependant. The variables appear more than once in the formal expression of the function f^. The CPU time for this 
inversion is about 255 seconds on a simple 1.67 Ghz Intel processor for a spatial surface resolution of 10^"*. 
In figure [8] one computes for this resolution the number of solution intervals according to candidates ones and the 
total surfaces ratio of the same intervals. The strong correlation between the two curves caracterizes the fact that the 
measures computations are well-defined because the probabilities are very close to the frequencies computed with ^p- 
algorithm independantly of the initial adjusments set. This confirms that the free algebra based arithmetic framework 
is relevant even for a semi-group IM calculations example. Due to the bissection, the algorithm computational 
complexity is exponential according to the iterations N, and it is not improved compared to SIVIA one. In our 
scheme, computational time is defined as 

Tcornp = 0(Af) = fc • 2^. (26) 

However, if the native function is differentiable on M", it is possible to define an adaptative mesh, with bissection 
spanned only in the space directions where the derivative magnitude is larger than a certain fixed value, because it 
is not useful to bissect in flat directions. This will obviously decrease the computational complexity constant k. 



V. CONCLUSION 



We have presented in this article an heuristic way to transfer real functions to inclusion ones depending on the space 
needed, semi-group or vector space of intervals. This permits to define a simple but very effic ient a lgorithm for set 
inversion inspired from SIVIA'^^, the ^-algorithm , which uses a free algebra based arithmeticP^'^si ^nd probability 
calculations. The convergence of this algorithm is assured for linear and non-linear problems. It offers several 
possibilities of applications, such as solving algebraic equations, differential equations, probability law of random 
variables calculations (discrete or continuous), topological analysis, numerical Lebesgue integrals computations, data 
analysis such as principal components analysis, and parameters identification. 
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